Testing Substellar Models with Dynamical Mass Measurements 



Trent J. Dupuy, a Michael C. Liu, and Michael J. Ireland 



1 Institute for Astronomy, University of Hawai'i, 2680 Woodlawn Drive, Honolulu, HI 96822 

2 School of Physics, University of Sydney, NSW 2006, Australia 



Abstract. We have been using Keck laser guide star adaptive optics to monitor the orbits of ultracool bina- 
ries, providing dynamical masses at lower luminosities and temperatures than previously available and enabling 
strong tests of theoretical models. We have identified three specific problems with theory: (1) We find that model 
color-magnitude diagrams cannot be reliably used to infer masses as they do not accurately reproduce the colors 
of ultracool dwarfs of known mass. (2) Effective temperatures inferred from evolutionary model radii are typi- 
cally inconsistent with temperatures derived from fitting atmospheric models to observed spectra by 100-300 K. 
(3) For the only known pair of field brown dwarfs with a precise mass (3%) and age determination («25%), the 
measured luminosities are ~2-3x higher than predicted by model cooling rates (i.e., masses inferred from Lboi 
and age are 20-30% larger than measured). To make progress in understanding the observed discrepancies, more 
mass measurements spanning a wide range of luminosity, temperature, and age are needed, along with more 
accurate age determinations (e.g., via asteroseismology) for primary stars with brown dwarf binary companions. 
Also, resolved optical and infrared spectroscopy are needed to measure lithium depletion and to characterize the 
atmospheres of binary components in order to better assess model deficiencies. 



1 Introduction 

Detailed theoretical models of stars, developed and obser- 
vationally tested over the last century, now underlie most 
of modern astronomy. However, the lowest mass stars (M* < 
O.1M ) are sufficiently cool (T eB < 3000 K) that the stan- 
dard, well-tested stellar models are not appropriate (e.g., 
due to dust formation in the photosphere). Objects below 
the hydrogen-fusing mass limit can cool to even lower tem- 
peratures as they have no sustained source of internal en- 
ergy generation. At temperatures below ~2000 K, dust and 
H2O are the major sources of opacity in the photosphere, 
and below ~1400 K methane absorption becomes impor- 
tant. Over the last decade, new theory has been developed 
to describe such low-temperature objects, which encom- 
pass brown dwarfs and gas-giant planets (e.g., JT|, fl2), 0, 
El, @, lfl3l ). These models now form the basis of our un- 
derstanding of all low-mass gaseous objects, from stars at 
the bottom of the main sequence to extrasolar giant plan- 
ets. Thus, rigorously testing them is vitally important. 

Dynamical masses for "ultracool" visual binaries (i.e., 
those with spectral types later than ~M7) are central to 
this effort, but such measurements have previously been 
impeded by observational limitations. Most ultracool bi- 
naries were discovered ~5 to 10 years ago, and the short- 
est estimated orbital periods are long (> 10-30 years). The 
corresponding physical separations (~l-3 AU) result in 
very small angular separations (< 0.2 arcsec) that can only 
be resolved using the Hubble Space Telescope (HST) or 
ground-based adaptive optics (AO). Direct distance mea- 
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surements via trigonometric parallax are also essential as 
they are needed to convert the observed angular scale of the 
orbit into physical units and equally importantly to provide 
a direct measurement of the luminosity. As a result, previ- 
ous to our work only three ultracool binaries had measured 
dynamical masses ([3|, [8|, |20|), and these were all rela- 
tively warm (> 2100 K). 



2 Dynamical Masses of Ultracool Binaries 

We have been using Keck laser guide star (LGS) AO di- 
rect imaging and aperture masking to monitor the orbits 
of ultracool binaries, enabling dynamical mass measure- 
ments for the lowest mass (30-75 Mj up ), lowest tempera- 
ture (1000-2800 K), lowest luminosity (10 -4 to 10~ 5 L ) 
objects known. This has allowed models to be tested in 
the unexplored area of parameter space shared by brown 
dwarfs and extrasolar giant planets. Keck LGS AO de- 
livers diffraction-limited imaging in the infrared for tar- 
gets over most of the sky (~0.05" FWHM; ~4x sharper 
than HST), so it is ideally suited for resolving short-period 
ultracool binaries. By performing a detailed analysis of 
these high-resolution images and accounting for small as- 
trometric shifts due to differential atmospheric refraction, 
we routinely achieve sub-milliarcsecond astrometric accu- 
racy and ~200 ^as for our best data. Such high quality as- 
trometry has allowed us to precisely measure binary orbits, 
with the error in the resulting masses (typically 3-10%) 
dominated by the uncertainty in the distance (Figure[TJ. 

We have also undertaken a substantial amount of sup- 
plementary analysis to enable these mass measurements: 
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Fig. 1. Our Keck LGS AO data combined with discovery and 
archival data from 5 to 10 years ago enables precise orbit deter- 
minations for ultracool binaries. Top left: HD 130948B and C are 
companions to a young solar analog (G2V, 0.8±0.2 Gyr), making 
them the first field brown dwarfs (L4+L4) with a well-determined 
age and masses (M tot = 0.109 ± 0.003 M Q ; (To)). Top right: 
2MASS J1534-2952AB is the first T dwarf binary with a dynam- 
ical mass (T5+T5.5, M tot = 0.056 ± 0.003 M G ), revealing incon- 
sistencies between the atmospheric model-derived temperatures, 
evolutionary model H-R diagram, and measured mass 1 14 1. Bot- 
tom left: LHS 2397aAB (M8+L7, M tot = 0.146±0.014 M Q ) is the 
first dynamical mass benchmark at the L/T transition, showing 
consistency between temperatures estimated from atmospheric 
and evolutionary models and supporting the idea that the temper- 
ature of the L/T transition is surface gravity dependent 1 1 1 1. Bot- 
tom right: 2MASS J2206-2047AB (M8+M8) is a pair of stars 
at the bottom of the main sequence that have ./-band colors 0.2- 
0.3 mag redder than predicted by evolutionary model tracks for 
objects of their measured masses |9|. This implies that masses 
and/or ages inferred from model color-magnitude diagrams will 
be in error for such objects. 



(1) We have re-analyzed archival HST images from 5 to 
10 years ago, improving astrometric errors by a factor of 
2 to 8 compared to published values, and this has proved 
critical for accurate orbit determination (e.g., IfTD . Ifl4l ). 

(2) We have developed a novel Monte Carlo technique to 
determine the orbital period probability distribution from 
motion observed between discovery and our first Keck data 
obtained > 5-10 years later (Figure |2). This has enabled 
us to accurately gauge target priorities and thus measure 
dynamical masses faster with a limited amount of tele- 
scope time. The orbital period, and thus the monitoring 
priority, of a visual binary is very uncertain from a sin- 
gle observation. In order to estimate the period probability 
distribution from a single observation, one must assume 
both a total mass and eccentricity probability distribution, 
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Fig. 2. We have developed a novel Monte Carlo technique to de- 
termine the orbital period probability distribution from orbital 
motion observed between only two epochs. Left: For the ultra- 
cool binary 2MASS J1728+3948AB, the relative positions of the 
A and B components are shown at the discovery epoch and ~6 
years later from our Keck LGS AO program (filled symbols). Us- 
ing the approach described in the text, we randomly drew orbits 
that pass through the two observed positions at the appropriate 
epochs (multi-colored lines). Right: The orbital period distribu- 
tion of the randomly drawn orbits (red) compared to the estimated 
period distribution using only the discovery epoch (black), fol- 
lowing the method of Torres 1191 , with ±lcr confidence limits 
shaded gray. Although the two-epoch distribution appears at face 
value to be broader and thus less precise, it is actually strongly 
preferred as it is free of the somewhat arbitrary assumptions re- 
quired by the single-epoch estimate (i.e., a uniform eccentricity 
distribution between < e < 1 and a total mass of 0.135 M Q ). 
In the case of 2MASS I1728+3948AB, we found that the orbital 
period is on the longer side (P > 26 years, 68.3% c.l.) of the 
original estimate (10^13 years, 68.3% c.l.), reducing its priority 
in our orbital monitoring program. 



and even with these assumptions the +lcr confidence lim- 
its span a factor of ~4 in orbital period jl9l . Our method 
utilizes the two positions and two times of two observa- 
tions taken several years apart to eliminate 4 of the 7 or- 
bital parameters (the two "geometrical" parameters semi- 
major axis and eccentricity; and the two "time" parameters 
period and time of periastron passage). This leaves only 
3 parameters (inclination, argument of periastron, and po- 
sition angle of the ascending node), which are just view- 
ing angles that we conservatively assume to be distributed 
randomly according to appropriate distributions. We use a 
Monte Carlo approach that results in an ensemble of or- 
bits that pass through the two observed positions at the ob- 
served epochs. The distribution of periods of these orbits 
is the period probability distribution. This method can, but 
does not necessarily, result in a narrower range of orbital 
periods; however, this method always results in a more ac- 
curate estimate of the period because it does not require an 
assumption about the eccentricity or total mass. 

(3) Finally, and most importantly, we have been conducting 
an infrared parallax program at the Canada France Hawaii 
Telescope (CFHT; Figure[3]). Precisely measured distances 
are critical for dynamical masses from visual binaries given 
the strong dependence of M tot oc cP . Only about 1 in 4 of 
the shortest period ultracool binaries have previously pub- 
lished parallaxes, so our program targeting ~30 binaries 
enables a greatly expanded sample of masses. 
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Fig. 3. Analysis of data from our infrared parallax program at 
CFHT demonstrates the capability to measure precise distances 
to visual binaries in our Keck LGS AO sample, which is criti- 
cal as the derived dynamical mass depends strongly on the dis- 
tance (Af tot °c d 3 ). Our parallaxes of "control" objects that have 
published parallaxes, such as the object shown here, agree well 
with published results. Our astrometric measurements are shown 
as open circles (error bars are smaller than the plotting symbols) 
along with the best-fit astrometric solution (solid black) which 
includes parallax (red dashed) and proper motion (blue dotted) 
components. 
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Fig. 4. We have tested substellar cooling rates by measuring the 
masses and luminosities of HD 130948B and C, brown dwarfs 
with a well-determined age of 0.8±0.2 Gyr from the primary star 
HD 130948A 1101 . Surprisingly, they are ~2-3x more luminous 
than evolutionary models predict. Two independent sets of the- 
oretical luminosity tracks are shown as colored isomass lines, 
where the thickness corresponds to the uncertainty in the mea- 
sured mass. Such a systematic error would imply that model- 
inferred masses are over estimated by ~ 20-30% in the usual 
case where only the luminosity and age are known (e.g., for di- 
rectly imaged extrasolar planets and determinations of cluster 
initial mass functions). The other possibility is that the age of 
HD 130948 A is inaccurately estimated despite using the latest 
gyrochronology calibrations 1151 . This can be resolved by ob- 
taining a better age estimate (e.g., from asteroseismology) and 
measuring masses and ages for more such systems. 



3 Testing Substellar Models 

Our Keck program has been yielding dynamical masses 
with the needed precision (3-10%) to perform fundamen- 
tal tests of theory (Figure [TJ. To date, we have identified 
three specific problems with substellar models: 

( 1 ) We find inconsistencies between predicted near-infrared 
colors and those observed for field objects of known mass 
over a wide range of spectral types ((9), iTTOl . IfTTl . |14|). 
For example, the M8+M8 components of 2MASS J2206- 
2047 AB and L4+L4 components of HD 130948BC appear 
to be aQ. 2-0.4 mag redder than model tracks, while the 
T5+T5.5 components of 2MASS J1534-2952AB appear 
to be »0. 2-0.4 mag bluer. This is likely due to imperfect 
modeling of dust for the late-M and L dwarfs and incom- 
plete methane line lists for the T dwarfs. Thus, masses 
and/or ages inferred from theoretical color-magnitude di- 
agrams should be treated with caution. 

(2) If theoretically predicted radii are correct, we have found 
that temperatures derived from atmospheric models are sys- 
tematically in error ((9), ifTUll . fl4'l). This points either to 
inaccurate theoretical radii or to incomplete modeling of 
such low-temperature atmospheres. The one object for which 
this is not the case is LHS 2397aB, currently the only mass 
benchmark at the L/T transition IfTTl . This is surprising as 
existing atmospheric models should not be appropriate for 
such transitional objects as their predictions are only valid 
in the limiting cases of maximal dust (Dusty, (6l) and the 
total absence of dust (COND, ||2)). Thus, in the case of 
LHS 2397aB it is likely that large systematic errors cancel 
out to produce apparent agreement. 



(3) For the only system with both a known mass and age, 
we have found the measured luminosities of the individ- 
ual components to be ~ 2-3 x higher than predicted (Fig- 
ure|3 1 10 1). This would imply that model-derived masses 
are significantly over estimated by ~20-30%. For example, 
an error of this magnitude is claimed to be needed to make 
the directly imaged extrasolar planets around HR 8799 dy- 
namically stable as their model-derived masses seem to be 
~30% too high fl2l . However, we emphasize that this ap- 
parent over-luminosity is based on a single system whose 
age is estimated from the primary star HD 130948 A. While 
this young solar analog is fortuitously amenable to multi- 
ple age-dating techniques using stellar rotation, chromo- 
spheric and x-ray activity, and isochrone fitting, precise 
age estimation is notoriously difficult for an arbitrary field 
star. Thus, refining the age estimate for HD 130948A (e.g., 
via asteroseismology) is essential. Also, more dynamical 
masses for the substellar binaries in triple systems with 
stars of known age (e.g., elnd Bab, G1417BC, GJ 1001BC) 
are critically needed to address this potential problem with 
model cooling rates. 

4 Future Work 

To make progress in understanding the problems we have 
found with substellar models, we are in the process of de- 
veloping a larger sample of masses spanning a wide range 
in temperature, mass, and age. For example, more mass 
measurements should determine whether the observed over- 
luminosity of brown dwarfs (FigureHJl persists for different 
surface temperatures. If so, this would point to a problem 
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with the interior structure model (e.g., convection) rather 
than a surface effect (e.g., magnetic fields). 

Resolved spectroscopy is also needed to enable detailed 
characterization of the individual components' tempera- 
tures, surface gravities, and abundances. Binaries with known 
masses can provide stronger tests of atmospheric models 
than field brown dwarfs of unknown age or mass. Such 
measurements will be able to assess poorly understood at- 
mospheric effects such as dust formation and sedimenta- 
tion (e.g., 0]) and vertical mixing that can drive nonequi- 
librium chemistry (e.g., Ifl6l ). 

In addition, optical spectra from HST/STIS will enable 
measurements of lithium absorption at 6708 A for bina- 
ries that we have shown are very close to the theoreti- 
cally predicted lithium burning limit at ^60 Mj up , provid- 
ing a novel test of substellar interior models ( IfTOl . IfTD ). 
As shown in Figure [5] independent groups make different 
predictions for the mass limit of lithium depletion. Bina- 
ries close to this limit offer the chance to empirically de- 
termine the lithium boundary for the first time, with a pre- 
cision comparable to our dynamical mass measurements 
(~3%). The theoretically predicted lithium boundary has 
never been directly tested in such a manner, even though it 
has been widely used, for example, to determine the canon- 
ical ages for young clusters such as the Pleiades |[T8l . 

Finally, a large sample of ultracool binary orbits will 
provide constraints on formation models for low-mass stars 
and brown dwarfs as the orbital parameters of binaries (e.g., 
semimajor axis and eccentricity) record the dynamical im- 
print of their formation process. For example, competing 
models predict very different eccentricity distributions for 
ultracool binaries: Stamatellos & Whitworth predict high 
eccentricities (<? > 0.5) for ultracool binaries formed via 
gravitational fragmentation in a massive circumstellar disk 
ITTTI , while Bate predicts modest eccentricities (e < 0.5) 
for ultracool binaries formed in a cluster environment Q . 
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